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ABSTRACT 

We present an observational study towards the young high-mass star forming region 
G23.44-0.18 using the Sub millimeter Array. Two massive, radio -quiet dusty cores MMl and 
MM2 are observed in 1.3 mm continuum emission and dense molecular gas tracers includ- 
ing thermal CH3OH, CH3CN, HNCO, SO, and OCS lines. The ^^CO (2-1) line reveals a strong 
bipolar outflow originated from MM2. The outflow consists of a low- velocity component 
with wide-angle quasi-parabolic shape and a more compact and collimated high- velocity 
component. The overall geometry resembles the outflow system observed in the low-mass 
protostar which has a jet-driven fast flow and entrained gas shell. The outflow has a dynam- 
ical age of 6 X 10^ years and a mass ejection rate --- 10~^ M© year~^ A prominent shock emis- 
sion in the outflow is observed in SO and OCS, and also detected in CH3OH and HNCO. We 
investigated the chemistry of MMl, MM2 and the shocked region. The dense core MM2 
have molecular abundances of 3 to 4 times higher than those in MMl. The abundance ex- 
cess, we suggest, can be a net effect of the stellar evolution and embedded shocks in MM2 
that calls for further inspection. 
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1 INTRODUCTION 



^7j The outflows take place in high-mass star forming regions at very 
^ early evolutionary stages (e.g. Birkmann et al.|20()6| [Beuther &[ 



|Sridharan|2007|[Longmore et al.|201l| . They inject large amount 
of hot gas and kinetic energy into the molecular cloud, causing 
intense shock waves that dramatically alter the chemistry of the 
surrounding environment {van Dishoeck & Bl ake 1998). How- 
ever, the morphological and dynamical properties of the out- 
flows on smaller scales, as well as their chemical effect to the 
young high-mass stellar cores are still to be further characterized. 



The high-mass star forming region G23.44-0.18 (023.44 
hereafter) was previously detected as a group of strong CH3OH 
masers i Walsh et al.| 19 981 which indicate the presence of young 
massive stars. It has a trigonometric-parallax distance of 5.88 
kpc (Brunthaler et al.|20'09) . In this letter, we report an observa- 
tion of the massive dusty cores MMl and MM2 associated with 
the CH3OH masers, as well as a high- velocity, intense bipolar 
outflow revealed in CO (2-1). The outflow is causing a shock 
emission and may be affecting the chemistry in MM2. Section 
3 presents the observational results. A further discussion on the 
outflow is given in Section 4. 1 and 4.2. The properties of the dusty 
cores are discussed in Section 4.3. A summary is given in Section 
5. 



2 OBSERVATIONS AND DATA REDUCTION 

The observational data is taken from the released SMA 
data archivq^ The observation toward 023.44 was car- 
ried out in September 2008. The phase tracking center is 
RA.(J2000)=18^34^39.25^ Dec.(J2000)=-8°3l'36.2^^ The point- 
ing accuracy of the SMA is ~ 3'^ The observation employed 8 
antennas in their compact configuration. The correlator has a 
total bandwidth of 4 OHz, centered at 220 OHz (LSB) and 230 
OHz (USB). The on-source integration time is 80 min. The SMA 
primary beam size (field of view) at this waveband (1.3 mm) is 
55'^ The correlator has a channel width of 0.812 MHz (1.1 km 
s~^). The system temperature is between 120 to 150 K. QSO 
3c454.3 and Neptune were taken as bandpass and flux calibra- 
tors, respectively QSO 1733-130 and 1911-201 were observed 
in every 30 and 50 minutes respectively during the observation 
sequence to track the antenna gains. The interferometer array 
has a shortest baseline of 17 m, corresponding to an u-v coverage 
for structures smaller than 20'^ 

The calibration and imaging were performed in Miriad^. 
The absolute flux level has an uncertainty of ~ 10%. The contin- 
uum map was extracted from the line-free channels in LSB. As 
the visibility data are INVERTed to the image domain, the syn- 
thesized beam size is 4.1^' x 3.7'^ 

Images of the 023.44 region at other wavelengths were also 
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Figure 1. The 1.3 mm continuum emission observed with the SMA 
overlayed on the IRAC composite image (blue=3.6, green=4.5, and 
red=8.0 p-m) image. The image center is RA.(J2000)=18^34"^39.25^ and 
Dec.(J2000)=-8°3l'36.2^'. The contour levels are -5, 5, 15, 25, 35, 45, and 
55 a (0.013 Jy beam~i). The synthesized beam (4.1" x 3.7") with PA = T 
from the north to NW. The squares denote the two strongest CH3OH 
masers. The crosses are the centers of the IRAC sources. 



Table 1. Physical parameters of the cores. 



parameter 


MMl 


MM2 


Unit 




m/^ 


12.315 ±0.049 


>15.8 


mag 




m^^ 


11.702 ±0.028 


>15.1 






m^,« 


11.435 ±0.023 


>14.3 






F(3.6M.m)^ 


1.42 ±0.28 


0.63 ±0.08 


mjy 




F(4.5^m)^ 


25.65 ±2.23 


14.48 ±1.96 






F(5.8M.m)^ 


47.84 ±2.22 


53.16±2.11 






F(8.0|am)^ 


32.92 ±2.04 


40.64 ±2.86 






F(450|Lim)^ 


27.3 ±1.5 


38.6 ±1.5 


Jy 




F(1.3mm) 


0.65 ±0.03 


1.42 ±0.03 






rrot(CH3CN) 


65 ±10 


110±10 


K 




rrot(CH30H) 


66±5 


60±6 






Mass 


120 ±30 


140 ±10 


Mo 




iV(H2) 


1.9 ±0.2 


2.5 ±0.2 


10^3 cm" 


-2 



Note. The flux densities at IR and sub-millimeter wavebands are 
taken from ^2MASS, ^Spitzer/IRAC, ^JCMT/ SCUBA. The remaining 
data all come from the SMA observation. 

extracted from several public archives. The Spitzer/IRAC im- 
ages and the point-source catalogue are taken the database of 
the GLIMPSE sky survey in the NASA/IPAC Infrared Science 
Archive^ We looked for the near-infrared counterparts in the 
2MASS point-source catalogue^. The continuum image at JCMT 
450 |j,m was also obtained to measure the flux densities. They are 
taken from the Canadian Astronomy Data Center (CADC) repos- 
itory of the SCUBA Legacy Fundamental Object Catalogu^ 



3 RESULTS 

Figure 1 shows the 1.3 mm continuum emission (white contours) 
overlayed on the RGB image of the three IRAC bands. It reveals 



two dust cores aligned from the north to south, denoted as MMl 
(northern core) and MM2 (southern core). Deconvoluted with 
the beam size, MMl has an mean diameter of 7.3^' (0.2 pc) while 
MM2 is more elliptical, with an extent of 9.7^' x 6.7^' (0.27 x 0.19 
pc) and position angle PA - 45° northeast. The strongest 6.7 GHz 
CH3OH masers ^Walsh et al.|1998) coincide with the centers of 
MMl and MM2. We referred to the MAGPIS 6 cm surve^Qfor 
the potential radio continuum emission from the ionized gas. It 
turns out that neither MMl nor MM2 has detectable emission 
above the sensitivity level of 2.9 mJy, indicating the both cores 
being prior to forming an Ultra-Compact Hii(UC Hii) region. The 
fluxes (or magnitudes) of MMl and MM2 at near- to far-Infrared 
are drawn from the 2MASS and IRAC point-source catalogues. 
The 450 |am fluxes are measured directly from the JCMT images. 
The obtained values are shown in Table 1. 

The 7 = 2-1 transition of ^^CO, ^^CO and C^^O and a num- 
ber of molecular species tracing dense gas are detected in our 
sidebands, including HNCO (10o,io-9o,9), SO (65-64), OCS (19- 
18), CH3CN \2Yi - lli^ (i^ = to 6), and four thermally excited 
CH3OH lines. In the online material, their spectra and physical 
parameters are provided in Figure SI and Table SI. The 2MASS 
Ks and JCMT 450 |j,m images are presented in Figure S2 and S3. 

Figure 2a shows the integrated images of the C^^O, CH3CN, 
HNCO, CH3OH (8o-7i), SO, and OCS lines. The C^^O has abroad 
distribution, with the emission peaks reasonably coincident with 
the dust cores. The other species show more compact emissions 
closely associated with the dust cores. The C^^O emission has 
comparable intensities for MMl and MM2, while the emissions 
of other dense-gas tracers are much brighter in MM2. Since the 
dust continuum and C^^O emission all closely follow the H2 dis- 
tributions, the intensity contrast suggests that MM2 have higher 
molecular abundances than MMl, as discussed in Section 4.3. 

A strong bipolar outflow is observed in ^^CO (2-1). Figure 
2b presents the CO (2-1) image integrated in the four velocity 
intervals. The outflow can be divided into a pair of low- velocity 
component (LVC) and high-velocity component (HVC) based on 
their distinct morphologies. The HVC roughly has \v\ ^ 20 km 
s~^ from the systemic velocity (101 km s~^). The bulk of the out- 
flow emission is in the southeast-northwest direction, with a po- 
sition angle of PA - 40° northwest. Figure 2c shows the beam- 
averaged SO and OCS spectra at MMl and MM2 center, and the 
shocked region (Aa = -ll'^ A(5 = 10^^. The ^^CO emission shows 
a similar feature with ^^CO despite being weak in the HVC. We 
present the channel maps of ^^CO and ^^CO in Figure S4 and S5, 
respectively. 



4 DISCUSSION 

4. 1 The physical properties of the outflow 

The presence of intense bipolar outflow in G23.44 suggests that 
the central stars are likely being formed via a disk-mediated ac- 
cretion ^ZhangetaL|2007| . The HVC is well collimated along the 
outflow axis, while the LVC, especially its red wing, has a more 
opened parabolic shape with a decreasing brightness from the 
edge to the center. The overall geometry of the outflow is remi- 
niscent of the outflow system in the low-mass protostar HH211 
(Gueth & Guilloteau|1999) . HH211 consists of a high-velocity jet- 
driven flow and a entrained low- velocity gas shell. The G23.44 
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Figure 2. (a) The velocity- integrated images of the molecular lines (contours) overlayed on the 1.3 mm continuum (gray). The integration is from 95 to 
110 km s~^ For C^^O (2-1), the contours are -4, 4, 8, 12... 40 a (1.35 Jy beam'^ km s'^). For CH3CN (122 - ll2)> the contours are -4, 4, 8... 24 a (0.6 Jy 
beam~i km s"^). For HNCO (10o,io -9o,9)> the contours are -4, 4, 6, 8... 20 a (0.34 Jy beam"^ km s"^). For CH3OH (80 -7i), the contours are -4, 4, 8... 48 a 
(0.5 Jy beam~i km s"^). For SO (65 - 54), the contours are -4, 4, 8, 12... 32 a (0.6 Jy beam"^ km s"^). For OCS (19-18), the contours are -4, 4, 5... 11 a (0.3 Jy 
beam~^ km s~^). The cross denotes the position of the shocked clump at (Aa = -11^', A(5 = 10^'). (b) The components of the ^^CO (2-1) outflow. The low 
velocity component (LVC) is integrated in (83,93) and (113,123) km s~^ for its blue- and redshifted lobes, respectively, while the high velocity component 
(HVC) is integrated in (65,75) and (140,170) km s~^ for the two lobes. For each lobe, the contours are -10, 10, 20, ..., 90% of the peak value. The dashed 
line denotes the axis of the outflow as indicated by the HVC. (c) The beam-averaged SO and OCS spectra at the center of MMl and MM2, and the shocked 
clump. The vertical lines denote the systemic velocity of i^isr = 101 km s~^ . 



Table 2. Abundance of the detected molecules. 



species 


A^x 


/x 


fa 
Jx,2ll 


(X) 


(cm-2) 






MMl(-1^^10") 






C^^O 


3.16±0.05(16) 


1.7±0.1(-7) 


- 


CH3OH 


2.5 ±1.2(15) 


1.3±0.6(-8) 


- 


CH3CN 


5.2 ±0.6(14) 


2.7±0.5(-9) 


- 


HNCO 


3.8 ±0.2(14) 


2.0±0.4(-9) 


- 


SO 


2.5 ±0.2(14) 


1.3±0.3(-9) 


- 


OCS 


3.7 ±0.2(14) 


1.9±0.3(-9) 


- 


MM2 a", 


-3") 






C^^O 


7.30 ±0.07(16) 


2.9±0.1(-7) 


1.7 ±0.2 


CH30H 


9.8 ±0.8(15) 


3.9±0.9(-8) 


3.0±1.5 


CH3CN 


1.6 ±0.2(15) 


6.3±0.7(-9) 


2.3 ±0.8 


HNCO 


2.20 ±0.05(15) 


8.8±1.5(-9) 


4.4 ±1.4 


SO 


1.45 ±0.06(15) 


5.8±0.9(-9) 


4.5 ±1.4 


OCS 


1.20 ±0.05(15) 


4.8±0.8(-9) 


2.5 ±0.7 


outflow-shock (-11", 10") 






C^^O 


2.4 ±0.07(16) 


= 2.0(-7)^ 


- 


CH3OH 


2.8 ±0.5(15) 


2.3±0.5(-8) 


- 


HNCO 


4.1 ±0.2(14) 


8.3±0.4(-10) 


- 


SO 


3.5 ±0.3(14) 


3.1±0.3(-9) 


- 


OCS 


8.0 ±0.3(14) 


6.6±0.3(-9) 


- 



a. The abundance ratio of MM2 to MMl. 

h. Assumed as the typical value in ISM. Abundances of the other 

molecules are then calculated by comparing to N{C^^O). 



outflow may represent a scaled-up version of this system, despite 
its collimation and shell structure being much less perfect than 
HH211. However, the size of the outflow in HH211 is only -10"^ 
AU, while the G23.44 outflow (10^ AU) has a 10-time larger spatial 
extent. On such a large scale, the outflow structure might begin 
to deteriorate due to increased turbulence in its environment. 
The CO outflow, especially the LVC also exhibits a weak 



emission feature in opposite direction to the major one, i.e. a 
blueshifted emission in the NW and a redshifted one in the SE. 
This can be a second flow by chance aligned in the same position 
angle with the major one. Alternatively, once the outflow direc- 
tion is close to the plane of the sky and its inclination angle be- 
comes smaller than the opening angle (~ 30°), its front and back 
sides will exhibit opposite Doppler shift. This scenario is well 
demonstrated by [Yen et aL](201 Figure 10 therein). We specu- 
late the second case to be more plausible. 

Assuming a local thermal equilibrium (LTE) and low opti- 
cal depth for the ^^CO line -wing emission, the CO column den- 
sity and the outflow mass can be inferred following the approach 
of |Garden et al.|(T991 1. The excitation temperature of the outflow 
is speculated to be slightly lower than in the dense core MM2 
(~ 70 K, Section 4.3). A value of Tgx = 50 K is adopted here. The 
outflow mass is estimated in each 1 km s~^ interval, then the 
total mass, momentum, and kinetic energy are calculated from 
M = Hyiniv), P = llym{v)v, and E^. = llym{v)v^/2. In calcula- 
tion, we assume an inclination angle of 6 = 20° to correct the 
velocity projection. The outflow turns out to have M ^ 5.5 M©, 



P ^ 360 Mq km s~^,Ej^^9x 10' 



46 



erg. 



Each of the blue and red wings in the HVC is resolved into 
two major clumps (Figure 2b, lower panel) denoted Bl, B2 and 
Rl, R2. These clumps in the outflow may indicate an episodic 
mass ejection (e.g. Q iu et al.|20 09j. The average central distance 
is 24^' (0.68 pc) for the pair of Bl-Rl, and 5.0^' (0.14 pc) for B2- 
R2. The two pairs both show a reasonable geometrical symmetry. 
However, the clumps are not exactly aligned across the MM2 cen- 
ter, instead with an offset to its west. This may indicate that the 
driven agent, or the stellar disk system is also displaced from the 
MM2 center. The dynamical age of the outflow t^y^ is estimated 
by dividing its spatial extent with the typical velocity {u = 40/ sinO 
km s~^). That yields r^yn = 6.1 x 10^ years for the Bl-Rl pair and 
f(jyn = 1-3 X 10^ years for B2-R2. Adopting the first value, we cal- 
culate the mass ejection rate to be M = ^/^dyn - l-^ ^ 10~^ Mq 
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yr ^ . One can see that the outflow is young and performing an 
intense mass ejection. 

4.2 Chemical enhancement in the outflow and shock 

In both low- and high-mass young stellar objects the out- 
flows have led to a rich chemistry in CH3OH, HNCO, and sul- 
phide s (Blake et al.||1987| [van der Tak et a'L]|2003| |j0rgensen| 



|et al.||2004| [Rodriguez-Fer nandez et al.||2010] etc.). In G23.44, 
the SO and OCS emission shows a noticeable emission clump 
at offset=(-ll'^10'0 (Figure 2a). The clump is almost perfectly 
aligned with the HVC axis of the CO (2-1), and reasonably co- 
incides with the peak R2, despite being closer to MM2 center. 
This strongly suggests that the clump traces a shock within the 
jet /outflow The clump is also weakly detected in HNCO and 
CH3OH (spectra shown in Figure SI). Their column densities 
are estimated using Equation (1) and (2) in Section 4.3. In cal- 
culation a temperature of Tgx = 50 K [same with that of CO (2- 
1)] is adopted, and the C^^O abundance is taken to be a con- 
stant of 2 X 10~^ as the typical ISM abundance which is close 
to /ci8o measured in MMl and MM2. The abundances of other 
molecules are then inferred from their intensity ratio relative to 
C^^O (2-1). The abundances of the shocked clump are listed in 
Table 2. 

In the low-mass protostars, the outflows often lead to a rich 
chemistry. In a few cases, the column density of HNCO and S- 
bearing molecules can been enhanced for orders of magnitude 
(e.g. LI 1 57 and L1448, Rodriguez-Fernandez et al.|2010|[lMfala et| 
al.|2Q10) . In comparison, the shocked region at (-ll'^ lO'O shows 



moderate abundances which are comparable or even slightly 
lower than in MM2. The enhancement level can be limited by the 
temperature and density in the shocked region, as well as the effi- 
ciency of the previous grain-surface chemistry. In addition, since 
the outflow has a short tdyn' the shock chemistry may still be un- 
dergoing and have not yet attained its maximum level. 

The blue lobe of CO (2 - 1) outflow are devoid of the SO and 
OCS emissions. As a possible explanation, the S-bearing species 
(e.g. H2S and CS as their progenitors) would be rather tightly 
bound to the dust grains (Blake et al.|1987) , and and a threshold 
shock velocity may be required for their desorption. Higher res- 
olution and primary shock products like SiO and H2S may better 
reveal the shock chemistry. 



4.3 The dense core properties 

The dust cores MMl and MM2 coincide with dark patches on the 
extended 8 ^im emission (Figure 1). The two cores are slightly 
connected with each other, suggesting that they might be the 
fragmentation products from the same natal cloud. 

The CH3CN (12j^ - 11^) lines provide a temperature diag- 
nosis for the cores (Figure 3a). We take an LTE radiative trans- 
fer model with a uniform rotation temperature and gas den- 
sity (Chen et al.|2006| [Zhang et al.|2007) to reproduce the ob- 
served CH3CN spectra at the dust-core center. We found the best 
fit at Trot = 65 K for MMl and Trot = 110 K for MM2. The mod- 
eled spectra ofK = and 1 towards MM2 is much lower than the 
observed line profile while the K > 1 lines are all well fitted. The 
excess emission in K = and 1 lines may largely arise from the 
cold outer-envelope gas. 

The rotation temperature is also evaluated from the ther- 
mally excited CH3OH lines at the core center using the rotation 
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Figure 3. (a) The solid line is the CH3CN (12^^ - 11^^) towards the cen- 
ter of MMl (top panel) and MM2 (bottom panel). The dashed line is the 
spectrum of the best fit from the LTE gas model, (b) The rotation diagram 
of the CH3OH lines. The filled circles are for the observed transitions. The 
vertical bars indicate the 3a errors. The linear least-squares fit is shown 
as the solid line. 



diagram method (Goldsmith & Langer|1999l|Liu et al.|2001||Biss-| 
chop et al.|2007) . The column density of the upper level A^u is 
related to the line intensity in the form 

Nu 3/c r 

gu 87r3vSju2 J ^ 

where gu is the upper-level degeneracy. S is the line strength, jj is 
the dipole moment. Nu depends on the total column density Nj 
in the form 

^u 
^ (2) 

where the partition function Qrot can be adopted from |Blake et| 
[aL]|1987| . Trot is then estimated from a least-square fit (Figure 
3b). The derived values are listed in Table 1. 

The Trot derived from CH3OH and CH3CN are remarkably 
similar for MMl. For MM2, however, the Trot of CH3OH is lower 
than that of CH3CN. As a probable reason, the CH3OH lines gen- 
erally have lower E^ than CH3CN (12^^^ - 11^^), then the cooler 
outer-part gas would contribute more significantly to the CH3OH 
lines than to the CH3CN (in particular K = 2to6).ln addition the 
CH3OH can be easily released from the dust grains by the shocks 
while the CH3CN production requires subsequent gas-phase re- 
actions and stellar heating (van Dis hoeck & Blake|1998j . CH3CN 
may therefore be distributed closer to the MM2 center (as seen 
in Figure 2a), and showing a higher Trot- 

For the other species, an independent temperature estima- 
tion is currently unavailable. We adopt Tqx = 60 K for MMl and 
Tex = 70 K for MM2 assuming an LTE excitation, the column den- 
sity of the molecules are calculated from their spectral lines to- 
wards the dust-core center, using equation (1) and (2). 

Assuming optically thin 1.3 mm dust continuum emis- 
sion, the total dust- and- gas mass can be estimated using Sy = 
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Mcore'KvBviT^)lgD'^. In the formula Sy is the total flux den- 
sity at 1.3 mm. g = 100 is the dust-gas mass ratio. ByiT^i} = 
{2hv^/c^)/[exp{To/T^] - 1] is the Planck function (Tq = hv/k). 
The dust opacity is set to be Xy = ^:230GHz = 0-9 cm^ g~^ of MRN 
dust grains with thin ice mantles at rin^ 



Ossenkopf& 



|Henning|1994 j. In the expression of Sy, using the beam-averaged 
flux density towards the core center, and dividing the obtained 
mass with rrtu^ , one can also calculate the peak column density 
N(H2) at MMl and MM2. The derived values are presented in 
Table 1. MM2 has a comparable but slightly higher N(H2) than 
MMl. 

The molecular abundances are estimated using fx = 
N{X)INiU2)' The /x and its ratio between MM2 and MMl are 
listed in Table 2. For the high- density tracers, the abundances in 
MM2 are generally 3 to 4 times higher than those in MMl. The 
SO shows the highest abundance ratio of 4.5. We note that the 
temperature adopted for MM2 (70 K) is conservative. A higher 
Tex would yield a lower N(H2) and higher N{X) for most species, 
hence leading to even higher abundances. For instance, adopt- 
ing the value of TrotCCHsCN) = 1 10 K as the Tqx for MM2, we will 
have N(HNCO) = 6.7 x 10^^ cm^ that is 3 times larger than the 
current value. 

The molecular distribution appears to be influenced by the 
outflow. Besides the coherent elongation of SO and OCS with the 
HVC of ^^CO emission (Figure 2b), the CH3OH and HNCO emis- 
sions are also largely distributed to the west of MM2. In addition, 
all the molecular lines have broader line widths in MM2 than in 
MMl, e.g. The SO line has Ai/pwHM = 11-3 in MM2 and 3.4 km 
s~^ in MMl. The line widths of MM2 can be broadened by a ki- 
netic energy input from the outflow. Considering its influence to 
the molecular line profiles, the outflow may also contribute to 
the higher abundances in MM2. Specifically, there would be un- 
derlying shocks which are induced as the outflow is interacting 
with the envelope and/ or infalling gas jArce & Sargent|2004l|Wu| 
|etal.|2009| . 

At the near- to mid-IR and sub -millimeter wavebands, the 
two cores have similar flux densities (Table 1). At IRAC bands 
MM2 is brighter than MMl at 5.8 and 8 ^tm but weaker at 3.6 and 
4.5 \im, while at 2MASS /, H, Ks bands, only MMl is significantly 
detected (Figure S2). Having a spectral energy distribution (SED) 
obviously skewed to the shorter wavelengths (A < 4.5 |j,m), MMl 
may have a more evolved stage for its stellar disk system ^Ro^ 
Ibitaille et al.|2006t . However, since the outflow in MM2 is close 
to the plane of the sky, it is also possible that a disk/toroid struc- 
ture lies edge-on in MM2 thus would largely absorb the short- 
wave emissions. A higher resolution will be helpful in revealing 
the morphology of the stellar system and underlying shocks in 
MM2. 



5 SUMMARY 

We present a multi-wavelength study towards the young high- 
mass star forming region G23.44-0.18. Two dusty cores, MMl 
and MM2 are observed in 1.3 mm continuum emission and the 
molecular lines of CH3OH, CH3CN, HNCO, SO, and OCS. The 
both cores have a pre-UC Hii evolutionary stage. A strong bipo- 
lar outflow is arising from MM2 as revealed by CO (2-1). The out- 
flow consists of a collimated, bi-polar high-velocity component 
and a more extended, parabolic flow at lower velocities. A clump 
of shocked gas is observed in SO (65 - 54) and OCS (19 - 18) lines 
and also detected in HNCO and CH3OH. The outflow shows a 



high momentum and mass-loss rate, and a short timescale of 
^dyn ~ 6 X 10^ years. 

In MM2 the high-density tracers CH3OH, CH3CN, HNCO, 
SO, and OCS show abundances of 3 to 4 times higher than those 
in MMl. To explain this difference, we suggest a combined effect 
of the stellar heating and underlying shock interactions in MM2. 
A further study with higher excited lines and improved resolu- 
tion can be performed to probe the chemistry, possible disk-jet 
system and mass accretion close to the MM2 center. 
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Table SI. Parameters of the molecular lines. 



Frequency 


Species 




^Isr 


^peak 


A^FWHM 


(GHz) 


(X) 


Transition 


(kms~i) 


(K) 


(kms-1) 


MMl (-1", 


10") 










219.5603... 


ci^o 


2-1 


102.2 


11±0.4 


3.2 ±0.5 


219.7982... 


HNCO 


10o,lO-9o,9 


101.4 


1.8±0.1 


3.6 ±0.4 


219.9494... 


SO 


65-54 


101.3 


4.2 ±0.1 


3.4 ±0.4 


220.0784... 


CH3OH 


80,8 -7i,6 


101.5 


1.0 + 0.3 


4.2 ±0.6 


220.7302... 


CH3CN 


I22-II2 


101.6 


1.3 + 0.2 


5.5 ±0.8 


229.7588... 


CH3OH 


8-1-70 


101.8 


3.9 ±0.1 


3.5 ±0.5 


230.0270... 


CH3OH 


3-2,2-4-1,4 


102.1 


1.3 ±0.3 


2.9 ±0.6 


231.0609... 


OCS 


19-18 


102.1 


1.9 ±0.2 


4.5 ±0.4 


231.2810... 


CH3OH 


102,9-93,6 


101.3 


0.6 ±0.2 


2.3 ±0.4 


MM2(l",- 


3") 










219.5603... 


ci^o 


2-1 
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231.0609... 
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6.7 ±0.4 


231.2810... 


CH3OH 


102,9-93,6 


102.1 


1.4 ±0.3 


4.3 ±0.3 


Shock (-11 


",10") 
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Figure SI. The molecular lines detected in the observational sidebands at the three positions. 
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Q 



20 



blue: 2MASS 2.2 micron 
red: IRAC 4.5 micron ^. 

white contours: SIVIA 1 .3 mm 
square: 6.7 GHz methanol masers 



20 -20 

RA offset (arcsec; J2000) 

Figure S2. The near to mid-IR RGB-coded image of G23.44. The blue is the 2MASS Ks, the green and red colors are IRAC 3.6 and 4.5 micron, respectively. 
The image center is RA.(J2000)=18^34"^39.25^ and Dec.(J2000)=-8°3l'36.2^'. The contour levels are -5, 5, 15, 25, 35, 45, and 55 a (0.013 Jy beam~i). The 
synthesized beam (4.1" x 3.7") with PA = 7° to the west. The squares are the two strongest CH3OH masers. The crosses are the centers of the IRAC sources. 




20 -20 

RA offset (arcsec; J2000) 



Figure S3. The SCUBA 450 \im image of G23.44. Contours are 10 to 90 per cent of the peak brightness (34.5 Jy beam ^). The beam size is 7.5". The cross 
denotes the center position RA.(J2000)=18^34"^39.25^ and Dec.(J2000)=-8°3l'36.2". 
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Figure S4. The ^^CO (2-1) channel images. The contours are -10, 20 to 90 percent of the peak intensity, which is 45 Jy beam ^ for the velocity range of 
[63,90] km s-^ and 18.9 Jy beam-^ for v= [110, 152] km s~^. Gray is the SMA 1.3 mm continuum emission. 
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Figure S5. The ^^CO (2-1) channel images. The contours are -10, 20 to 90 percent of the peak intensity which is 12.2 Jy beam~^ for the velocity range of 
[84,93] kms~i and 8.8 Jy beam" ^ for v= [108, 117] kms~^ 



